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ABSTRACT 

A self-similar solution for time evolution of quasi-spherical, self-gravitating accretion flows is obtained 
under the assumption that the generated heat by viscosity is retained in the flow. The solutions are 
parameterized by the ratio of the mass of the accreting gas to the central object mass and the viscosity 
coefficient. While the density and the pressure are obtained simply by solving a set of ordinary differential 
equations, the radial and the rotational velocities are presented analytically. Profiles of the density and 
the rotational velocities show two distinct features. Low density outer accreting flow with relatively 
flat rotation velocity, surrounds an inner high density region. In the inner part, the rotational velocity 
increases from the center to a transition radius where separates the inner and outer portions. We show 
that the behaviour of the solutions in the inner region depends on the ratio of the heat capacities, 7, 
and the viscosity coefficient, a. 

Subject headings: accretion, accretion discs - hydrodynamics 



1. INTRODUCTION 

Accretion processes are now believed to play a major 
role in many astrophysical objects, from protostars to 
disks around compact stars and AGN. Such systems have 
been studied at different levels depending on their physical 
properties. Geometry of the disk (thin or thick), transport 
of the thermal energy inside the disk, self-gravity of the 
accreting gas and the magnetic fields are among the most 
important factors which shape any theory for such sys- 
tems. For simplicity, traditional models of accretion disks 
assume geometrically thin configuration and neglect self- 
gravity of the accreting material (Shakura and Sunyaev 
1973; Pringle 1981). These models have been extended by 
considering large scale magnetic fields, polytropic equation 
of state and better understanding of the mechanism of an- 
gular momentum transport (e.g., Livio and Pringle 1992; 
Ogilvie 1997; Duschl, Strittmatter and Biermann 2000). 
However, the key ingredient in all such models is that the 
generated heat by turbulent viscosity does not remain in 
the flow. In other words, all of the viscously dissipated 
energy was assumed to be radiated away immediately. 

Another type of accretion flow known as Advection- 
Dominated Accretion Flow (ADAF) has been proposed, in 
which generated heat by viscosity can not escape from the 
system and retains in the flow (Ichimaru 1977; Narayan 
and Yi 1994). During recent years, ADAFs have been 
paid attention as plausible states of accretion flows around 
black holes, active galactic nuclei or dim galactic nuclei (for 
review see Kato, Fukue and Mineshige 1998). Using sim- 
ilarity method, Ogilvie (1999) (hereafter; OG) extended 
original steady state self-similar solutions to time depen- 
dent case. His solutions describe quasi-spherical time de- 
pendent advection dominated accretion flows. 

On the other hand, many authors tried to study the 
effects related to the disk self-gravity. Recent observa- 
tions show significant deviations from Keplerian rotation 
in objects that are believed to be accretion disk, and in 
some cases, there are strong evidences that the amount of 



mass in the disk is large (e.g., Drimmel 1996; Greenhill 
1996). The role of self-gravity in accretion disks was dis- 
cussed by Paczyhski (1978), who studied vertical structure 
of the disk under the influence of self-gravity. Some au- 
thors investigated role of self-gravity on waves in the disk 
(Lin and Pringle 1987; Lin and Pringle 1990). Within the 
framework of geometrically thin configuration, Mineshige 
and Umemura (1996) extended classical self-similar ADAF 
solution (Narayan and Yi 1994) and found a global one di- 
mensional disk solutions influenced by self-gravity both in 
the radial and perpendicular directions to the disk. In 
another study, Mineshige and Umemura (1997) extended 
previous steady solutions to the time dependent case while 
the effect of self-gravity of the disk had been taken into ac- 
count. They used isothermal equation of state and so their 
solutions describe a viscous accretion disk in slow accretion 
limit. Then, Mineshige, Nakayma and Umemura (1997) 
extended this study by obtaining solutions for polytropic 
viscous accretion disks. Also, Tsuribe (1999) studied self- 
similar collapse of an isothermal viscous accretion disk. 

Although solution of Mineshige and Umemura (1996) 
describes self-gravitating ADAFs, applicability of the so- 
lution restricts to geometrically thin configurations. The 
most appropriate geometrical configuration for advective 
flows, rather than thin, is quasi-spherical as has been con- 
firmed by many authors (e.g, Narayan and Yi 1995). OG 
presented the first semi-analytical quasi-spherical, time- 
dependent ADAFs solution. Here we note that although 
OG considered advection-dominated accretion flow in a 
point mass potential, it is straightforward to relax this as- 
sumption so as to describe the effect of self-gravitation on 
quasi-spherical advection-dominated accretion flows using 
time-dependent similarity solutions. 

This paper is organized as follows. In section 2 the gen- 
eral problem of constructing a model for quasi-spherical, 
self-gravitating accretion flow is defined. The self-similar 
solutions are presented in section 3, and the effects of the 
input parameters are examined. We summarize the results 
in section 4. 
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2. FORMULATION OF THE PROBLEM 



3. SELF-SIMILAR SOLUTIONS 



We start with the approach adopted by Ogilvie (1999) 
who studied quasi-spherical accretion flow without self- 
gravity. In his approach, the equations are written in 
spherical coordinates (r, 6, ip) by considering the equato- 
rial plane 9 = j and neglecting terms with any 9 and ip 
dependence. It implies that each physical variable repre- 
sents approximately spherically averaged quantity. So, all 
the physical quantities depend only on the spherical ra- 
dius r and time f. Self-gravity of the accreting gas can 
be considered simply by Poisson equation and the corre- 
sponding term in the momentum equation. The governing 
equations are the continuity Equation, 

t + ^>'> = °- m 

the Equations of motion, 

^L +v ^L + l^P + ^l + ^Ml = V l ( 2 ) 
dt r dr p dr dr r 2 r ' 

- ^ll^l&l (3) 



y dt dr r 2 dr 

the Poisson's Equation, 
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and the energy Equation, 
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2r 9 



(5) 

In order to solve the equations, we need to assign the 
kinematic coefficient of viscosity v. Although there arc 
many uncertainties about the exact form of viscosity, au- 
thors introduce some prescriptions for v regarding dimen- 
sional analysis or just based on phenomenological consider- 
ations (e.g, Duschl, Strittmatter and Biermann 2000). We 
employ the usual a prescription for the viscosity, which we 
write as (Shakura and Sunyaev 1973) 

(6) 



v — a 



pfl k ' 



pp,p^pp, V r><p 



where fik = (GMj/r 3 ) 1 ' 2 is the Kcplerian angular velocity 
at radius r. This prescription originally introduced for vis- 
cosity in thin accretion disk, however, it has been widely 
used for studying dynamics of thick accretion disk, such 
as ADAFs. 

To simplify the equations, we make the following sub- 
stitutions: 

vv rv ,^ — > **,r — > rr,t — > it, 

(7) 
where 

GM+ f . 1 „ - . GM* ,^ 

r t v47tG/j r 
Under these transformation, the continuity equation does 
not change and the rest of equations are cast into 
dv T 
~dt 



dv r 
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dr p dr dr 
, d d ad 
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r 2 dr dr r ' 
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r 2 dr dr 
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(11) 
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We look for self-similar solutions of Equations (1) and 
(9)-(12) and reduce this system into a set of ordinary dif- 
ferential equations. We define a self-similar variable 



(t -t) 2 / 3 ' 

where t < to and demand that 



(13) 

(14) 
(15) 
(16) 
(17) 
(18) 



Substituting these expressions into the above Equations, 
we obtain the following set of self-similar Equations: 
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(19) 
(20) 
(21) 
(22) 



Interestingly, the continuity Equation is integrable and 
gives ^ 2 R(^+2>V/2) = C, where C is a constant and should 
be zero. Thus, 



(23) 



Also, we can obtain rotational similarity function $(£) 
simply by integrating Equation (22), 



where 



3 7 



3a 1-7 



(24) 



(25) 



and $ s is rotational velocity at some £ s . For simplicity, 
we assume that £ s defines the outer boundary of the ac- 
cretiong gas. This Equation correctly gives $(£ = 0) = 0, 
and reaches to a maximum $ m by increasing £ from zero 
to a point at £ m . One can simply show 
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where £ = 3$ S /4A and 



4/3 



A 



a/3 



(26) 



(27) 



Clearly, behaviour of rotational velocity and the other 
physical variables except for the radial velocity, depend on 
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£ s and $ s . We can consider both £ s and <& s as free param- 
eters, however, it's possible to determine them uniquely, 
if we parameterize self-similar solutions using conserved 
quantities, e.g. mass of the system. As can be seen from 
the definition of A, the above self-similar solutions are ap- 
plicable within the range of 1 < 7 < 4/3 if we consider a 
positive value for a. 

We can derive asymptotic solutions when approaching 
the origin £ = as follows 



By inserting dS/d£, from Equation (19) into the above 
Equation, we finally obtain 

9 P 2 dP 
M = M^l + g£ s 3 - f s (^) s 1]. (41) 

Another relation is obtained from Equation (31), as follows 
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(28) 
(29) 
(30) 



These asymptotic solutions are valuable when performing 
numerical integrations to obtain similarity solutions start- 
ing from x — > + . Note that the behaviour of the solutions 
near to origin only depends on A and a. 

We now proceed to solve the rest of above similarity 
Equations. By substituting Equation (24) into Equation 
(20), we obtain 



4 q d£_ 



where 



/(0 



|^[3^ + 4A(V/4_ 4 -3/4 0] _ 
6a\ £ s I 



(31) 



(32) 



After some algebra, from this Equation and Equations (19) 
and (21), the following Equation is obatined: 

15 „„d 2 P 

(33) 

where 
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dt 



s 3 = 
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(34) 
(35) 
(36) 
(37) 



This Equation can easily be integrated considering ap- 
propriate boundary conditions. We can parameterize the 
solutions by £ s as a function of the mass of the accret- 
ing gas and the viscosity parameter a. For simplicity, we 
assume P(£ s ) = i? s and P(£ s ) = P B . Total mass of the 
accreting gas is 



M = 4ir pr dr. 
Jo 

Using similarity solutions, this Equation reads 



= M* f 
Jo 



RCd£_. 



Substituting from Poisson Equation (21), we have 



(38) 



(39) 



(40) 
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/y,g(3$ s -2A^). 



(42) 
(43) 



Now, we can use the standard fourth-order Runge-Kutta 
scheme to integrate nonlinear ordinary differential equa- 
tion (33) from sufficiently small £ to £ s . Given £ s and 
$ s and using asymptotic solutions near to the origin (i.e., 
equations (29) and (30)), one can start numerical integra- 
tion and equation (41) gives corresponding ratio of masses, 
i.e. M/M*. Before presenting results of integration, it 
would be interesting to explore the typical behaviour of 
rotational velocity $(£). 

Just as an illustrative example we follow another ap- 
proach: lets assume that we know i? s , P s and M/M* from 
numerical integration and so, one can determine £ s and $ s 
from Equations (41) and (42) uniquely. For example, if we 
set P s — and (dP/d^) s — and consider non-zero but 
negligible P s , the outer boundary of the system is deter- 
mined analytically 



e. = ^2A* + 4.5(^ + 1) -^^ + 4.5(^ + 1)), 

(44) 

(45) 



and 



In this case, Figure 1 shows £ s and $ s as functions of a 
for 7 = 5/4 and M/M+ = 0.6, 1, 1.5 and 2.5. As this Fig- 
ure shows for fixed ratio of masses, the outer radius of 
the system increases by increasing the viscosity coefficient 
a. Also, if the ratio of masses increases, the outer radius 
increases irrespective of the exact value of a. However, 
rotational velocity at £ s (i.e. <J> S ) is not very sensitive to 
values of a or M/M* for large values of a. If the viscosity 
coefficient a tends to small values, on the other hand, the 
rotational velocity $ s increases significantly. 

In Figure 2, we plot rotational velocity of representa- 
tive cases with a = 0.01,0.1 and M/M* = 0.6,1.5,2.5. 
As this Figure shows, rotation of the flow increases from 
the center to the outer radius. There are two regimes as 
rotation of the flow concerns. While in the outer part, pro- 
file of the rotational velocity in similarity space is nearly 
flat, in the inner portion this velocity strongly increases. 
However, for low values of a (e.g. 0.01), rotational velocity 
reaches to a maximum then with nearly constant and small 
slope decreases. Also, the rotational velocity increases as 
the parameter a decreases and the outer region with flat 
rotation profile becomes larger, as the ratio of the mass 
of the accreting gas to the mass of the central object in- 
creases. In other words, while rotational behaviour of the 
inner part is nearly independent of the ratio of masses, the 
outer portion's rotation and it's extension are sensitive to 
that ratio. 



4 



In Figure 3, typical behaviour of the density and the 
pressure in similarity space are shown. We can see in 
the inner region the solutions can clearly be described by 
asymptotic solutions (29) and (30), however, outer region 
has different density and pressure profiles. Interestingly, 
behaviour of solutions in the inner region is independent 
of the extension of the system (£ s ) and the rotational ve- 
locity $ s . As the ratio of M/M* increases, the extension 
of the outer region becomes larger. Flows with large val- 
ues of a have lower central mass concentration comparing 
to accretions with low values of a. Generally, one can 
say quasi-spherical, self-gravitating accretion flows of our 
model consist of two parts, one inner portion with high 
density and an outer part with nearly lower density and 
flat rotation profile. However, in both regions the radial 
velocity is the same, irrespective of the input parameters. 

From Equations (16) and (17), we can simply show that 
the ratio of the radial velocity to the rotational velocity 
is v r (r,t)/v v (r,t) = V/$>. Figure 4 shows this ratio for 
a = 0.1 and 0.01 and three different values of M/M*. 
Clearly, this ratio of velocities are much smaller than unity, 
implying that at each radius r and at any time t, rotational 
velocity is greater than the radial velocity. However, as 
the viscosity parameter a increases, the ratio increases as 
well. This behaviour is easy to understand; because we 
showed that radial similarity velocity V is proportional to 
similarity variable £, independent of the input parameters. 
However, rotational similarity velocity $ decreases by in- 
creasing the parameter a, as can be seen from Figure 2. 
So, the ratio V/<& increases in the case that a increases. 

As we showed the total mass of the accreting gas is 
conserved. We can simply write the integral represent- 
ing the angular momentum as J = 4ir J pv v r^dr. Using 
the similarity solutions, we can show that the total angular 
momentum J is proportional to (to — t) 1 ^. In non-self- 
gravitating quasi-spherical accretion flow, J is conserved 
and the total mass M is proportional to (to — t)^ 1 ^ 3 as 
has been shown by OG. But our solutions imply that the 
total mass of accreting material is conserved and angular 
momentum is decreasing. Also, the central density p c in- 
creases as p c cx (to — t)~ 2 . For non-sclf-gravitating flow, 
we see another behaviour: p c cx (to — i)~ 7 / 3 . However, the 
radius of the flow in both cases decreases in proportion to 
(h-tf/\ 



4. SUMMARY 

In this paper, we have studied quasi-spherical accretion 
flow, in which heat generated by viscosity retained in the 
flow. In opposition to the usual studies performed up to 
now, we have considered self-gravity of the flow. We de- 
rived similarity solutions for such flows which are appli- 
cable within the range < 7 < 4/3, if we consider posi- 
tive value of a. Radial and rotational velocities have been 
obtained analytically. Obtained solutions parameterized 
by the ratio of the disk mass to the central object mass, 
M/M+, and the viscosity parameter, a. We showed that 
the extension of the accreting gas depends on this ratio 
and the viscosity parameter. More importantly, these in- 
put parameters have direct effect on the rotational veloc- 
ity. 

Our solutions are different from the solutions by OG 
in various respects. We found that the radial similarity 
velocity is in proportion to £, implying no critical point. 
This fortunate circumstance let us to integrate rest of the 
Equations simply, although physically one should bear in 
mind this kind of velocity (i.e., independent of input pa- 
rameters) is as a result of mathematical limitations of sim- 
ilarity method. At the outer edge of the accreting gas, the 
density and the pressure have low values comparing to the 
central region. 

Other viscosity laws have been proposed, for instance 
the /3-prescription which is based on analogy with tur- 
bulence observed in laboratory sheared flows and gives 
v cx v v r (Duschl, Strittmatter and Biermann 2000). Al- 
though we have not explored self-gravitating accretion 
with this prescription, there are self-similar solutions with 
similarity indices the same as those have been found with 
a-prescription in this paper. However, ordinary differen- 
tial equations governing the similarity physical variables 
are different and should be solved numerically. It would 
be interesting one compare the similarity solutions with 
/^-prescription with those we have obtained here. 



I thank the anonymous referee for his/her report which 
showed an error in the initial version of the paper, prompt- 
ing me to find the correct form of the asymptotic solutions 
near to the origin. 
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Fig. 2. — Self-similar rotational velocity corresponding to 7 = 5/4 and R s = P E = 0.01. Each curve is marked by the ratio of the mass of 
accreting gas to the mass of the central object, M/M+, and £ s and <3> s as (M/M+, £ s , <E> a ). 
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Fig. 3. — Self-similar density and the pressure corresponding to 7 = 5/4 for different values of viscosity coefficient, i.e. a = 0.01 and 0.04. 
Curves are labeled the same as Figure 2. 




Fig. 4. — Ratio of the radial velocity to the rotational velocity for 7 = 5/4 and R B = P B = 0.01 and a = 0.1, 0.01. Curves are labeled by 
the ratio of the masses, i.e. M/M+ 



